^^AD-A0B2  130 

TELEDYNE  BROWN  ENGINEERING  HUNTSVILLE  AL 
THERMAL  RESPONSE  MODELS  MBALL» (U) 

OPTICAL  SYS— ETC  F/«  20/13 

SEP  79  WE  SMITH 

DASG60-79-C-0050 

UNCLASSIFIED 

MN-S0-9-79-893 

NL 

MH-SB-9-TS-H3 


THERMAL  RESPONSE  MODEL:  MBALL 


September  1979 


DT1C 

jELECTEJ 
k  MARI  21980. 


I 


DISTRIBUTION  STATEMENT  * 

Approved  for  public  r«lea»e| 
Distribution  Unlimited 


WTELEDYNE 
BROWN  ENGINEERING 

Cummings  Research  Park*  Huntsville.  Alabama  35807 


IS  lfO 


LETTER  OF  TRANSMITTAL  MSE79-BMDATC-4735 


^TELEDYNE 

BROWN  ENGINEERING 


TO: 


FROM: 


SUBJECT: 

DATE: 


Director 

Ballistic  Missile  Defense 
Advanced  Technology  Center 
F.  0.  Box  1500 
Huntsville,  Alabama  35807 

Attn:  Max  Hardwick,  ATC-D 

Optical  Systems  Department 
Systems  Division 
Teledyne  Brown  Engineering 
Cummings  Research  Park 
Huntsville,  Alabama  35807 

Transmittal  of  MB ALL 

7  November  1979 


To  enable  more  faithful  signature  calculations  of  decoy  target 
concepts  applicable  to  the  Optical  Discrimination  Program,  TBE  has 
developed  a  modification  (MBALL)  to  the  Optical  Signatures  Code  (0SC) . 
This  modification  supercedes  all  previous  versions. 

The  following  is  a  list  of  materials  to  be  distributed  as  speci 
fled  on  the  attached  page: 


•  MBALL  manual 

•  OSC  Listing  (including  MBALL  in  EXOHEAT) 

•  A  short  description  of  the  new  6-degree-of- 
freedom  modification  to  the  trajectory  program 
BALL IS 

•  Update  decks  of  MBALL  and  the  6-degree-of- 
freedom  option  for  the  OSC  VI  cycle. 


If  you  have  any  questions  concerning 
signed  at  205-532-1355. 


Warren  Smith 


the  above,  contact  the  under- 


1 


LETTER  OF  TRANSMITTAL  MSE7 9-BMDATC-4 735 


Director 

Ballistic  Missile  Defense 
Advanced  Technology  Center 
P.  0.  Box  1500 
Huntsville,  Alabama  35807 

Attn:  Max  Hardwick,  ATC-D 

Ballistic  Missile  Defense  System 
Command 

P.  0.  Box  1500 
Huntsville,  Alabama  35807 

Attn:  BMDSC-C 

Defense  Technical  Information 
Center 

Cameron  Station 
Alexandria,  Virginia  22314 

McDonnell  Douglas  Corporation 
5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Attn:  Herdman 

H.  Ginell 

MIT  Lincoln  Laboratory 
P.  0.  Box  73 

Lexington,  Massachusetts  02173 

Attn:  Wade  Kornegay 

Nichols  Research  Corporation 
71910  South  Memorial  Parkway 
Suite  A 

Huntsville,  Alabama  35802 

Attn:  C.  Horgen 

Sandia  Laboratories 
P.  0.  Box  5800 

Albuquerque,  New  Mexico  87185 
Attn:  Bruce  Bulmer 


Enclosures: 

Letter  of  Transmittal 
Manual 


Letter  of  Transmittal 


Letter  of  Transmittal 
Manual  (2) 


Letter  of  Transmittal 


Letter  of  Transmittal 

Manual 

Listing 

Card  Decks 


Letter  of  Transmittal 

Manual 

Listing 

Card  Decks 


Letter  of  Transmittal 

Manual 

Listing 

Card  Decks 


2 


BALLIS  EXO  6 DOF  ANALYTIC  SOLUTION 


Insert  into  OSC  VII  Basic  Manual: 

BALLIS  6DOF  ADDITION:  IKIND  -  4 

In  order  to  save  computer  time  for  exoatmospheric  6  degree  of  freedom 
trajectories  where  two  moments  of  inertia  are  equal,  an  analytic  solution 
to  Euler's  equations  with  no  external  torques  has  been  developed.  Such 
a  solution  is  valid  for  axi-symmetric  bodies  above  the  top  of  BALLIS 's 
atmosphere  (•vlO6  feet) .  The  rationale  for  this  modification  is  to 
reduce  the  computer  time  required  by  BALLIS.  The  same  input  is  required 
for  this  option  as  for  option:  IKIND  “3,  IKIND6  “  0.  The  option  is 
accessed  by  inputting  IKIND  as  4  on  Card  4A.  Cards  4B,  5  and  6  are 
also  required. 

Insert  into  Table  4-11  Basic  Manual: 


IKIND 


■4  Exo-axisynmetric  6D0F  Analytic  Solution 
Cards  5  and  6  required. 


/  ij  MN-SD-9- 79-893 

J  ^ 


THERMAL  RESPONSE  MODEL: 


MBALL  / 
. ^ 


J' 


Sept  «rtw -3:979 


| 


Sponsored  By 

BALLISTIC  MISSILE  DEFENSE  ADVANCED  TECHNOLOGY  CENTER 
DEPARTMENT  OF  THE  ARMY 
HUNTSVILLE,  ALABAMA 


Prepared  By 


OPTICAL  SYSTEMS  DEPARTMENT 
SYSTEMS  DIVISION  ' 
TELEDYNE  BROWN  ENGINEERING 
HUNTSVILLE,  ALABAMA 


l)  o  i  y  /7 


ABSTRACT 


MBALL  is  a  subroutine  of  EXOHEAT  for  surface  temperature  calcu¬ 
lations  for  convex  objects  of  light  to  intermediate  thermal  mass  (in-depth 
heat  conduction  is  not  important).  It  includes  surface  conduction  and 
internal  radiative  coupling  between  stations,  using  the  natural  fluxes 
(Sun,  albedo,  earthshine)  calculated  in  EXOHEAT.  MBALL  must  be  used 
with  the  OSC  VI  BASIC  option. 


APPROVED  BY: 


J.VV.  Beaupre 
Manager 

Optical  Signature  Analysis 


Section 


N.  A.  Passino 

Deputy  Manager 

Optical  Systems  Department 


TABLE  OF  CONTENTS 


1. 

2. 

3. 


4. 


5. 

6. 
7. 


Page 


INTRODUCTION  .  1-1 

1.1  MBALL  Overview  . .  1-1 

1.2  MBALL  Capabilities .  1-1 

FUNCTIONAL  ANALYSIS .  2-1 

THEORY .  3-1 

3.1  Heat  Equation  .  3-1 

3.2  Solution  of  the  Heat  Equation  .  3-21 

3. 3  Open  Surface  Option  .  3-24 

INPUT  SPECIFICATIONS .  4-1 

4.1  Disc  Utilization .  4-1 

4.2  Card  Input .  4-1 

4.3  Accessing  MBALL  in  the  BASIC  Option  .  4-9 

EXAMPLES .  5-1 

5.1  Hater-Jacketed  Balloon . . .  5-1 

5.2  Replica  . 5-13 

HINTS  AND  DIAGNOSTICS  .  6-1 

REFERENCES  .  7-1 


I 


i 


J 


LIST  OF  ILLUSTRATIONS 


Figure  Title  Page 

1- 1  OSC  VI  Program  Structure .  1-2 

2- 1  Exoheat  Program  Structure  .  2-2 

3- 1  Target,  Solar,  Earth  Geometry  .  3-6 

3-2  Nosecap  Geometry  .  3-8 

3-3  Station  Normals  (Exterior)  .  3-9 

3-4  Target  Geometry .  3-10 

3-5  Nosecap  Polar  Angles  . 3-11 

3-6  Radiation  Exchange  Geometry  .  3-12 

3-7  Internal  Radiation  Geometry  .  3-15 

3- 8  Conduction  Geometry  .  3-17 

4- 1  Specific  Heat  and  Enthalpy  of  Water  .  4-8 

5- 1  Balloon  Example  Input  .  5-2 

5-2  Balloon  Example  Output  .  5-6 

5-3  Replica  Geometry  .  5-14 

5-4  Replica  Example  Input  .  5-15 

5-5  Replica  Example  Output  .  5-16 


iv 


LIST  OF  TABLES 


Table  Title  Page 

1-1  MBALL  Capabilities  .  1-3 

A  A 

3- 1  Attitude ^Averages  for  Incidence  Factors:  G(n,s)  «* 

n*s  0(n*s)  •  . .  3-7 

4- 1  EXOHEAT  (with  MBALL)  Disc  Files  .  4-1 

4-2  RVSNTH  Input  Summary  .  . .  4-2 

4-3  EXOHEAT  (with  MBALL)  Card  Input  Summary  .  4-4 

4-4  OSC  Thermophysical  Property  Code  Words  .  4-9 

4- 5  Modeling  the  Heat  of  Fusion  of  Water  4-9 

5- 1  p,  0^,  k  for  Water  and  Teflon  .  5-5 


v 


.....  a  ..  -  - - »,.|i|(-  ■■ 


1.  INTRODUCTION 


1.1  MBALL  OVERVIEW 

Subroutine  MBALL  was  developed  to  provide  target  surface  tempera¬ 
tures  in  an  exoatmospheric  environment  for  balloon  and  replica  shapes  of 
light  to  intermediate  thermal  mass.  These  temperatures  are  used  by  the 
Optical  Signatures  Code  (OSC)  (Ref.  1)  to  generate  long  wavelength  infrared 
(LWIR)  signatures  for  threat  discrimination  analysis.  It  is  necessary 
that  vehicle  geometry,  deployment,  and  external  flux  information  be  sup¬ 
plied  by  the  OSC  VI  to  MBALL,  which  performs  temperature  calculations  for 
thermally  light  objects. 

The  OSC  execution  sequence  is  shown  in  Figure  1-1.  MBALL  is  located 
in  the  EXO/ENDO  thermal  response  block  as  a  subroutine  of  EXOHEAT. 

1.2  MBALL  CAPABILITIES 

Table  1-1  is  a  summary  of  MBALL' s  capabilities.  For  each  point 
of  the  vehicle  surface  at  which  temperatures  are  calculated,  MBALL  uses 
averages  of  in-depth  material  properties  (assuming  a  layered  skin  of  dif¬ 
ferent  materials)  and  computes  an  average  temperature  as  if  the  skin  were 
made  up  of  one  homogeneous  material.  This  temperature  is  then  assigned 
to  the  surface  of  the  vehicle  for  signature  generation. 
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FIGURE  1-1.  OSC  VI  PROGRAM  STRUCTURE 


TABLE  1-1.  MBALL  CAPABILITIES 


Targets 

•  Thin  to  Intermediate  Thermal  Mass 

•  Replica  or  Balloon  Shapes 

Thermal  Fluxes 

•  Solar 

•  Albedo  From  EXOHEAT  Subroutine 

•  Earthshine  (EXOENV) 

•  Molecular 

Station  Coupling 

•  Internal  Radiation 

•  Thermal  Conduction:  Longitudinal  and 

Transverse 

Thermal  Response 

•  Phase  Change  Capability 

•  Thermal  Properties  Updated  with  Temperature 

Data  Base 

•  Earthshine  Data:  Models  Based  on  NIMBUS 
Observations  (Ref.  2)/L0WTRAN4  Calculations 

•  Thermophysical  Properties:  OSC 

•  Thermophysical /Chemical  Data  Base:  OSC 

Options 

•  Radiative  Equilibrium  (Thermal  Mass  =  0) 

•  Open  Surfaces 

2.  FUNCTIONAL  ANALYSIS 


Target  surface  temperatures  are  determined  in  the  OSC  by  EXOHEAT 
for  exoatmospheric  conditions  (greater  than  400  kft)  or  by  RVTEMP  for 
endoatmospheric  conditions  (less  than  400  kft).  These  programs  need  tar¬ 
get  position,  velocity,  and  deployment  information  as  a  function  of  tra¬ 
jectory  time  to  compute  the  temperatures.  This  information  is  supplied 
by  BALLIS  or  by  the  user.  The  temperatures  are  then  used  by  BIDIREC  to 
generate  radiometric  signatures,  with  material  optical  properties  supplied 
by  SELECT.  BALLIS,  EXOHEAT,  RVTEMP,  SELECT,  and  BIDIREC  are  all  contained 
in  the  OSC  and  are  called  by  the  OSC  program  BASIC?,  depending  on  the  user 
input . 

MBALL  is  a  subroutine  of  EXOHEAT  and  replaces  EXOHEAT  in-depth 
temperature  calculations  (i.e.,  replaces  Subroutine  EXOTMP)  for  objects 
of  light  to  intermediate  thermal  mass.  Heat  flux  calculations  are  per¬ 
formed  in  EXOHEAT  (Subroutine  EXOENV)  for  a  given  ballistic  trajectory 
and  vehicle  geometry  to  determine  external  heating  rates  (from  Sun,  albedo, 
and  Earth  emission)  and  these  are  passed  to  MBALL,  bypassing  EXOHEAT ' s 
temperature  subprograms.  The  vehicle  geometry  required  by  MBALL  is  com¬ 
puted  in  RVSNTH,  a  subroutine  of  BASIC?.  Previously,  RVSNTH  was  used  only 
with  the  endoatmospheric  heating  routine,  RVTEMP,  but  has  been  modified 
to  supply  MBALL  with  the  necessary  parameters.  RVSNTH  does  not  calculate 
SELECT  or  BIDIREC  data  when  MBALL  is  called,  however,  so  the  inputs  to 
SELECT  and  BIDIREC  must  be  supplied  by  the  user.  The  calculations  that 
determine  the  external  heating  fluxes  are  described  in  the  EXOHEAT  manual 
(Ref.  3)  (written  before  the  MBALL  option  existed).  Reference  1  contains 
the  input  to  SELECT  and  BIDIREC.  Figure  2-1  shows  the  EXOHEAT  program 
flow.  A  description  of  how  MBALL  is  called  in  the  BASIC  option  is  in 
Section  4. 
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FIGURE  2-1 .  EXOHEAT  PROGRAM  STRUCTURE 
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3.  THEORY 


3.1  HEAT  EQUATION 

The  heat  equation  expressing  the  time  rate  of  change  of  tempera¬ 
ture  T  at  any  point  in  a  material  of  density  p,  thermal  conductivity  k, 
and  specific  heat  Cp,  with  an  internal  energy  source  density  S,  can  be 
written: 


p  Cp  |£  =  V  •  (k  VT)  +  S  (3-1) 

where  the  term  kVT  can  be  thought  of  as  an  energy  density,  whose  divergence 
gives  an  energy  source  density.  The  term  VT  is  the  spatial  gradient  of 
the  temperature. 

In  MBALL,  the  thermally  light  skin  of  a  replica  or  balloon  shape 
is  modeled  by  dividing  it  into  N  pieces  defined  by  the  user,  where  each 
piece  has  associated  with  it  an  outer  surface,  an  inner  surface,  and  up 
to  four  surfaces  that  adjoin  adjacent  pieces.  The  outer  surface  com¬ 
municates  with  the  external  environment  via  radiative  coupling,  the 
inner  surface  communicates  with  the  other  interior  surfaces  in  the  same 
manner,  and  the  adjoining  surfaces  communicate  through  heat  conduction. 

Each  piece  also  has  an  associated  density,  specific  heat,  and  thermal 
conductivity.  The  values  chosen  for  these  quantities  for  a  given  piece, 
say  the  jth  piece,  are  an  appropriate  average  of  the  exact  values  taken 
over  the  volume  of  the  piece.  By  integrating  Equation  3-1  over  the 
jth  piece  with  this  point  in  mind  (the  subscript  j  refers  to  average 
values),  one  obtains: 

f  Pj  cPj  sr  dvj  mf  -  °j  <*j 

y  Areas 

^  Bounding  Vj 

A 

where  nj  is  a  unit  normal  pointing  out  of  surface  dAj  of  piece  j,  and 
S  ■  0;  i.e.,  the  assumption  is  made  that  no  independent  sources  of  energy 
(batteries  and  wires)  exist  in  the  skin  of  our  model.  Because  the 
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material  quantities  are  averages,  they  can  be  taken  as  constant  over  the 
volume  of  Integration,  The  area  Integral  can  be  resolved  into  its  component 
parts : 


9Tj  f 

pj  cpj  vj  9 r  mJ  (kVT)  *  ^ 

Aj  outer  surface 


j  +  f  (k$T) 


nj  dAj 
Aj  inner  surface 


/  (k$T)  •  n^  dAj  . 
Aj  adjoining  surfaces 


(3-3) 


The  boundary  conditions  for  the  outside  and  inside  surfaces  of  the 
jth  piece  can  be  written  as  follows,  where  the  jth  piece  is  assumed  to 
radiate  as  a  greybody  with  temperature  Tj  and  emissivity  e  j : 


A 

ai  * 

k^T  - 

-£j 

+  Qext.j 

(outside) 

(3-4) 

*  - . 

k$T  - 

”ej 

°V  +  Qint,j 

(inside) 

(3-5) 

where  prime  marks  denote  internal  normals  and  emissivities ,  and  a  is  the 
Stefan-Boltzmann  constant.  Qext  j  *-s  the  average  power/area  incident 
on  the  exterior  of  piece  j  and  is  due  to  the  external  environment. 

Q^nt  j  is  the  average  power/area  striking  the  interior  of  piece  j  and 
is  due  to  the  other  radiating  pieces. 

Substituting  Equations  3-4  and  3-5  into  Equation  3-3  for  the 
appropriate  surface  integrals  yields: 


pj  cPj  VJ 


-/  (w,j  -£j  aV> 

Aj  outer 

+/  «w,i  -Ei  °v> 

Aj  inner 
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dA 


j  • 


A  j  :  adjoining 

(3-6) 
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Because  these  are  averages,  and  assuming  AQuter  “ 


(thin  skin) : 


3Ti 


pJ  CP1  VJ  ■  “W.j  *J  +  AJ  -  (ej  +  EJ>  01 . *  *J 


+  2  Aadjolning  ke£f  TfdJ°lnln6 

_ i  _ • _ _  iiA  ,  i  ,  _  ,  _  •  _  _ 


aj  oining 
areas 


j : adjoining 


(3-7) 


where  the  conduction  integral  of  Equation  3-6  has  been  approximated  by 
a  sum  over  all  adjoining  surfaces  of  the  piece  j.  The  terms  Aa<jj0ining* 

keff’  Adjoining*  and  ^.-adjoining  refer  t0  the  contact  area  between 
piece  j  and  its  neighbor,  a  properly  averaged  thermal  conductivity 

between  piece  j  and  its  neighbor  (discussed  below),  the  average  tempera¬ 
ture  of  the  adjoining  piece,  and  the  distance  between  the  centers  of  j 
and  its  neighbor,  respectively. 

It  is  convenient  to  rewrite  Equation  3-7  in  the  following  way 
to  simplify  the  bookkeeping: 

3T1  — 
pj  'pj  VJ  aT  *  W.j  Aj  +  Mij  Ti4+SKlJ  Ti<3'8) 

i  i 

where  the  matrices  M^j  and  K^j  contain  the  radiation  and  conduction 
terms,  respectively,  which  are  discussed  in  detail  below.  The  sums 
are  over  all  pieces,  including  the  jth. 

Because  of  the  need  to  model  materials  that  have  rapidly  changing 
specific  heats  as  a  function  of  temperature  (i.e.,  water  near  0°C),  it 
is  necessary  to  generalize  the  left-hand  side  of  Equation  3-8: 

9Tj  s  r  TJ 

pj  cpj  vj  St - -  pj  vj  -h  J  °vi  (Tj">  dTj~  <3-9> 
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where  TQj  Is  the  initial  temperature  of  piece  j ,  and  Tj  is  the  tempera¬ 
ture  after  some  time  t.  Thus,  the  equation  solved  for  Tj  by  MBALL  can 
be  written: 


»j  VJ  -k  [  1  <TP  +2  M« Ti' 

•/  rp 


Oj 


♦  2T  11 


(3-10) 


The  linearization  of  the  left-  and  right  -hand  sides  necessary 
for  an  iterative  solution  for  Tj  is  described  in  detail  in  Section  3. 

3.1.1  External  Sources 

There  are  four  natural  external  sources  of  radiation  considered: 
sunshine,  albedo,  earthshine,  and  molecular  collisional  heating.  These 
parameters  are  calculated  in  EXOHEAT  and  are  passed  to  MBALL.  Reference 
3  contains  a  complete  discussion  of  these  quantities. 

e  Sunshine  -  A  solar  exoatmospheric  irradiance  of  1353 
W/m2  is  assumed.  The  Sun  position  is  input  in  one  of 
three  ways: 

A  Subsolar  point  input:  latitude  and  longitude 

A  Subsolar  point  calculated  from  month,  day,  GMT 

A  Subsolar  point  calculated  from  month,  day,  local 
(24-hour)  time  at  a  given  longitude 

e  Albedo  -  The  reflectivity  (a)  of  the  sunlit  Earth  is 
assumed  to  be  0.4.  The  albedo  radiance  of  a  sunlit 
element  of  the  Earth  is 


Na 


a  S  cos  6 
IT 


(3-11) 


where  S  is  the  solar  constant  and  6  is  the  angle  between 
the  Earth  normal  and  the  direction  to  the  Sun. 
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•  Earthshine  -  Six  Earth  radiance  maps  are  available 
to  EXOHEAT,  corresponding  to  day /night  conditions 

for  winter,  summer,  and  equinox.  Each  map  corresponds 
to  a  5-  by  5-deg  latitude,  longitude  grid.  The  seasons 
are  determined  by  the  month,  and  Sun  position  determines 
the  time  of  day.  Reference  2  contains  a  magnetic  tape 
appendix  containing  earthshine  data.  The  earthshine 
map  may  be  modified  with  NAMELIST /SOLAR  input  (Ref.  3, 

P.  6) 

•  Molecular  -  At  altitudes  below  106  ft,  a  molecular 
collisional  heating  flux  is  used  of  the  value 

Qalr  ’  1  pair  V’  (3*12> 

where  Pflir  is  the  air  density  and  V  is  the  velocity 
of  the  target. 

The  effect  that  the  above  contributions  have  on  the  energy  incident  on 
a  given  area  element  of  the  target  depends  on  the  geometry  of  the  Sun, 
Earth,  and  target  configuration,  and  where  the  area  element  is  located 
on  the  target.  Only  the  flux  component  normal  to  the  area  element 
is  considered.  Reference  3  contains  a  detailed  description  of  how  this 
problem  is  treated.  Figure  3-1  shows  the  geometry. 

EXOHEAT  has  the  option  of  calculating  three  different  flux 
averaging  methods  (Table  3-1) : 

•  Instantaneous  -  Energy  incident  on  a  given  target 
element  is  that  which  is  determined  by  the  instantaneous 
geometry  (i.e.,  position  of  the  Sun  and  Earth,  and 
location  of  the  element  on  the  target) 

•  Roll-Averaged  -  Energy  incident  on  a  given  target 
element  is  the  normalized  average  of  the  energy  the 
element  would  see  through  a  full  revolution  about 
the  target  longitudinal  axis;  therefore,  all  elements 
lying  in  a  ring  centered  on  the  body  axis  experience 
the  same  incident  energy. 

•  4ir  Average  -  Incident  energy  is  the  normalized  average 
of  the  energy  that  the  element  would  see  through  a 
"tumble"  over  4ir  steradians.  All  elements  on  the 
target  receive  the  same  incident  energy. 
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earth  geometry 


TABLE  3-1.  ATTITUDE  AVERAGES  FOR  INCIDENCE  FACTORS: 
G(n,s)  =  n*s  0(n* s ) 


MODE* 

<  G  > 

COMMENT 

1 

G 

Not  averaged.  Instantaneous 
heat  flux  calculation 

2 

1  r  2ir 

2?  /  d$G 

Roll  averaged 

0 

3 

f  G  =  i 

4tt  average:  fast  random 
tumbling 

♦Input  parameter 


MBALL  has  the  capability  of  treating  transverse  as  well  as 
longitudinal  heat  conduction  as  described  below,  so  mode  1,  instantaneous 
heating  rate,  is  necessary  to  make  use  of  this  ability. 

3.1.2  Internal  Sources,  Vehicle  Geometry,  Radiation  Matrix 

MBALL  considers  radiative  coupling  between  the  inner  surfaces 
of  a  replica  or  balloon  shape.  It  is  possible  to  model  a  completely 
closed  structure,  or  one  in  which  the  baseplate  has  been  removed 
(see  Section  3.3).  Mo  internal  source  of  energy,  such  as  batteries 
or  wires,  are  included. 


The  user  has  the  choice  of  modeling  either  a  balloon  (sphere) 
or  an  axially  symmetric  replica  that  can  have  a  spherical  nosecap, 
up  to  three  frusta,  and  a  flat  baseplate.  It  is  necessary  to  input  the 
length  of  the  nosecap  such  that  it  meets  the  first  frustum  tangentially 
(Figure  3-2).  The  surface  area  of  the  vehicle  can  be  divided  into  as 
many  as  50  pieces,  called  stations,  where  each  station  has  an  outer  and 
inner  surface  normal,  an  initial  temperature,  an  average  thickness. 


Q  =  CONE  ANGLE  OF  FIRST  FRUSTUM 
LI  =  R1  (1-sin  Q) 


FIGURE  3-2.  NOSECAP  GEOMETRY 


specific  heat,  thermal  conductivity,  solar  absorptance,  outside 
emissivity,  inside  emissivity,  and  density  defined  by  the  user.  All 
stations  on  the  nose  section  (nose,  frustum,  and  base  are  sections)  have 
the  same  areas,  which  are  calculated  by  the  program  RVSNTH.  The  orienta¬ 
tion  of  the  outer  normal  of  each  station  is  defined  by  the  angle  it  makes 
with  the  z  axis  (looking  toward  the  nose)  and  the  roll  angle  <f>  (Figure  3-3). 
The  running  length,  RL,  distance  from  the  nose  along  the  axis  z,  and  the 
distance  from  the  axis,  r,  of  the  center  of  the  station  are  calculated 
in  RVSNTH  (Figure  3-4) .  The  polar  and  azimuthal  angles  are  input  by 
the  user  to  EXOHEAT  to  determine  the  external  heating  rates.  The  polar 
angles  for  the  normals  on  the  frusta  and  base  are  chosen  to  agree  with 
their  respective  slopes,  but  the  choice  for  the  polar  angles  on  the 
spherical  nose  stations  should  be  done  such  that  all  nose  stations  have 
equal  areas  (Figure  3-5) .  This  ensures  a  consistent  treatment  of  external 
and  internal  radiation  in  EXOHEAT  and  MBALL.  The  roll  angle  $  is  measured 
counterclockwise  from  the  vehicle  x  axis,  looking  from  the  nose  to  the 
base  (Figure  3-3) . 
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THIS  EXAMPLE  HAS  7  LONGITUDINAL  BY  4  AZIMUTHAL  =  28  STATIONS  TOTAL 


NUMBERING  OF  STATIONS 
INCREASES  LONGITUDINALLY: 


4>i :  STATION  1-7 
<f>2  8-14 

<J>3  15-21 

$4  22-28 


y 


©=  POSITION  OF  SURFACE  NORMAL 
(BASE  NORMALS  ARE  NOT 
VISIBLE) 

Z  IS  ALONG  AXIS,  OUT  OF 
PAGE 


VEHICLE  AS  VIEWED  FROM  THE  FRONT 


FIGURE  3-3.  STATION  NORMALS  (EXTERIOR) 
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r*:  DISTANCE  FROM  AXIS  (ft) 

Zf:  DISTANCE  FROM  NOSE  (ft) 

RLT:  RUNNING  LENGTH  (ft) 

A*:  STATION  AREA  (ft2) 

t.:  THICKNESS  (THK)  (ft) 

♦CALCULATED  IN  RVSNTH 

NOTE:  B*  IS  NOT  PASSED  TO  EXOHEAT 
BY  RVSNTH,  AND  MUST  BE 
CALCULATED  BY  THE  USER  FOR 
EXOHEAT  INPUT,  SEE  FIGURE 
3-5. 


FIGURE  3-4.  TARGET  GEOMETRY 


IN  THIS  EXAMPLE, 
THERE  ARE  THREE 
NOSECAP  STATIONS 
OF  EQUAL  AREA. 


FOR  N  NOSECAP  STATIONS: 

»,  •  »■-  ['  -  St 


FOR  A 

WHERE 


SPHERE  (BALLOON),  LI  *  2  ■  R1  , 

*  cos*:  I1  * 

N  *  TOTAL  NUMBER  OF  LONGITUDINAL 
STATIONS 


FIGURE  3-5.  NOSECAP  POLAR  ANGLES 


The  stations  can  be  distributed  over  the  sections  as  desired, 
with  the  only  restriction  being  that  the  number  of  roll  angles  times  the 
number  of  polar  angles  be  less  than  or  equal  to  50.  Thus,  it  is  possible 
to  have  50  longitudinal  stations  along  a  single  <j>,  or  50  azimuthal  stations 
at  a  single  longitudinal  position. 

There  is  a  small  ambiguity  in  the  term  "station"  as  used  in  the 
RVSNTH  input.  There,  it  means  the  number  of  stations  on  a  given  section, 
assuming  only  a  single  roll  angle.  The  program  automatically  multiplies 
this  by  the  number  of  roll  angles  to  get  the  correct  number  of  stations 
on  the  section. 

The  balloon  is  handled  similarly  to  the  replica,  keeping  in  mind 
the  criteria  of  equal  areas  for  the  determination  of  the  polar  angles. 

Section  4  describes  the  input  to  RVSNTH  and  EXOHEAT. 

The  radiation  matrix,  M* j ,  will  now  be  derived.  The  power,  d4P, 
incident  on  an  area  dAj  from  a  blackbody  source  of  area  dA^  can  be 
written  as  (see  Figure  3-6): 


FIGURE  3-6.  RADIATION  EXCHANGE  GEOMETRY 


°  Ti 


cos  ou  cos  on 


a^,  ctj  -  the  angles  the  normals  to  dAi,  dAj  make  with  ri  j , 
the  vector  pointing  from  dA^  to  dAj 

£j  -  the  emissivity  of  the  piece  j 

Ti  -  the  temperature  of  the  source 

0  -  the  Stefan-Boltzmann  constant. 

The  source  surface  over  which  Equation  3-13  is  integrated  is  the 
interior  surface  of  the  replica  or  balloon  model,  which  has  been  divided 
into  N  stations  as  described  previously.  Because  the  temperature 
and  internal  emissivity,  e' ,  of  each  station  are  average  values  over  the 
station,  the  integral  of  Equation  3-13  over  the  entire  inner  surface 
(excluding  r^j  =  0)  can  be  written  as  a  sum  over  each  station  surface 
with  e' ,  T  outside  the  integral,  which  results  in: 


(3-14) 


.  BB 
where  g  , 


the  blackbody  radiation  exchange  matrix,  is  given  by 


(3-15) 


To  take  into  account  the  fact  that  the  sources  on  the  interior 


are  not  blackbodies,  one  replaces  g  in  Equation  3-14  by 

T  ,  BBt  .  BBt, 

g  ■  lg  (1  -  p  g  )  E  ] 


(3-16) 


where  g  is  the  true  radiation  exchange  matrix  and  superscript  T  stands 
for  "transposed".  The  reflectivity  p'  is  given  by  p'  -  1  -  z' . 


For  a  closed  surface,  one  has 


0  for  the 


Holes/open  surfaces  are  correctly  treated  by  setting  cl  m  pi  ■ 
surfaces  transparent  to  radiation  in  the  calculation  of  g.  Equation 
3-17  no  longer  applies  to  g  if  some  surfaces  are  transparent. 

The  total  power  incident  on  Aj  is  approximated  by 

Pj  *  £j  Aj  aE  gij  Ti\-  (3"18) 

i 

The  key  to  the  calculation  of  g  is,  by  Equation  3-16,  the 

T}T>  BB 

calculation  of  g  .  The  calculation  of  gj^  is  accomplished  by  the  subdivision 
of  Aj  into  r)  longitudinal  strips  with  center  positions  rm,  and  a  numerical 
integration  is  performed: 

A  11 

gji  =  —  cos  °  cos  a  *m) / I2  •  (3-19) 

m=l 


To  ensure  that  no  serious  errors  result  from  this  numerical 

BB 

procedure.  Equation  3-17  is  imposed  on  g  . 

Excluding  external  sources  for  the  moment,  the  total  radiation 
flux  into  Aj  is  the  amount  incident  from  the  interior  minus  the  internal 
and  external  emission  of  Aj : 


Pj  TOTAL 


•  AJ 
(3-20) 


which  can  be  written  as: 


Pj  TOTAL  *  23  ^  Ti 

i 


(3-21) 
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INSIDE  SURFACE 
OF  A  TYPICAL 
FRUSTUM 


i  (CUT  INTO 
n  STRIPS) 


FIGURE  3-7.  INTERNAL  RADIATION  GEOMETRY 


where 


Mij 


-o(e, 


+  ep 


5iJ 


♦•2 


(3-22) 


with  ^ij  “  j0 

The  term  the  "self  contribution",  can  be  understood  as  a 

correction  to  the  energy  radiated  away  internally  by  a  curved  piece. 
One  edge  of  the  piece  radiates  back  toward  the  opposite  edge,  thus 
reducing  the  total  outward  energy  flow. 


where  j±0,  j±4>  label  the  pieces  adjacent  to  the  jth  piece,  in  the  longi¬ 
tudinal  and  azimuthal  directions,  respectively  (see  Figure  3-8).  The 
remaining  terms  are  defined  as: 


where 


M  -  T  <*i  +  «'j+4>  -  \  wj  +  *i-t> 


[(“•146  +  *M 

/RLi  + 

(*■1  1 

\  2  / 

\  2  /. 

V  2  / 

ix3±e,j  ■  *  (la3±9  -  V 
“*W  •  ) 


(3-25b) 

(3-25c) 

(3-25d) 

(3-25e) 


kefr(n,D 


where  n,  £  *  j,  j±6,  j±4>  (3-25f) 


and  where  r,  <J>,  t,  KL  are  defined  in  Figure  3-4. 

By  defining  longitudinal  and  transverse  conduction  coefficients 
of  the  form: 


a 


j±e 


V±) 


keff(j*8,i) 

AXj±8,j 


(3-26a) 


6j±4> 


A(T±)  keffW,j) 

AX^,j 


(3-26b) 


Equation  3-24  can  be  written  as: 


£  ^adjoining  k.ff 


ad j oining 


j  .'adjoining 


“  °j-e  Tj-e  +  aj+e  Tj+e 

+  ^-4,  Tj— 4>  +  Bj+<|> 

-  Tj  (aj-e  +  aj+e  +  Bj-<j>  +  BjV 

(3-27) 
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Thus,  by  defining  a  conduction  matrix  of  the  form 


•l4« 

*  *7*4  4  4 

^7*4 


0 


(3-28) 

The  sum  (including  i  *  j)  £  k^  T^  will  contain  all  conduction  terms 
associated  with  the  jth  piece,  so  that  Equation  3-27  can  be  written  as: 


E 

adjoining 


Aad joining  ^eff 


adjoining  ; 
^ j : ad j  oining 


E  *ij 


(3-29) 


MBALL  updates  the  K^j  with  the  changing  temperatures  of  the  stations  through 
the  thermal  conductivity  ke^,  which  is,  in  general,  a  function  of  the 
temperature . 

Because  the  end  pieces,  i.e.,  the  first  nose  piece  or  the  last 
baseplate  piece,  have  only  one  longitudinal  conduction  contribution,  the 
appropriate  elements  of  the  matrix  K^j  are  set  to  zero. 

3.1.4  Thermal  Mass 

The  left  side  of  Equation  3-10  can  be  written  as: 


’j  VJ  £  fi  <T3<T">  iT"  -/^  "j 


(3-30) 


where  mj  ■  Pj  cpj  is  defined  as  the  thermal  mass  of  the  jth  piece, 
and  is,  in  general,  a  function  of  the  temperature  because  it  contains 
the  specific  heat.  If  the  vehicle  being  modeled  has  a  skin  structure, 
i.e.,  if  it  is  built  up  of  layers  of  materials  with  different  densities 
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and  thermal  properties,  an  appropriate  average  of  these  quantities  must 
be  made  over  all  of  the  layers,  because  only  a  lumped  thermal  mass  is 
treated  in  MBALL  for  each  station.  Equations  3-31  through  3-33  give 
the  average  density,  specific  heat,  and  thermal  conductivity  that  are 
necessary : 


E 


£  at.  k 

k  J’ 


c 


p,j 


£  p3>k  AXj>k  cpj,k 

E  p3.k  AXj,k 

k 


kj 


AX.  ,  k.  , 
J,k 


(3-31) 


(3-32) 


(3-33) 


where  AX^  p..  cp^  and  kj  k  are  the  thickness,  density,  specific 
heat,  and  thermal  conductivity,  respectively,  of  the  kth  layer  of  piece 
j.  Because  the  specific  heat  and  thermal  conductivity  are  input  by 
the  user  for  different  temperatures,  an  appropriate  cpj  and  k^  should 
be  input  for  each  temperature.  The  temperature  calculated  for  the  jth 
piece  by  MBALL  is  the  average  temperature  over  the  thickness  of  the  skin. 
Rapidly  varying  specific  heats  are  accounted  for  in  the  solution  of  the 
heat  equation  described  in  the  next  section. 

The  option  exists  to  set  the  thermal  mass  to  zero  '(BALL  ■  2) 
for  thin  skins,  in  which  case  the  temperature  of  each  station  is  solved 
for  assuming  radiative  equilibrium; 

(cj  +  epot,*  -  Qext>J  +  Qint>J  <«*> 


where  Qext  j*  Q^nt  j  are  the  external  and  internal  power  per  unit  area, 
respectively.  Incident  on  the  jth  piece. 
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where  At  is  the  time  step,  T0j  is  the  initial  temperature  of  piece  j,  and 
Tj'  is  an  estimated  value  for  the  new  temperature  Tj  after  the  time  At. 

The  determination  of  At  and  T"  is  discussed  below.  Because  the  temperature 
change  is  being  sought  over  a  finite  interval  At,  it  is  appropriate  to 
replace  the  instantaneous  value  of  the  right-hand  side  of  Equation  3-10 
by  a  value  averaged  over  the  interval: 


aj  +  £  (M«  V  *  Ti>  —  +  T>  Ai 

i 

+  2  (M*j  +  Kij  (3-36) 

i 

where  T^,  the  average  temperature,  is  given  by  »  1/2  (T0^  +  Tj) ,  with 
T0i»  Ti  the  old  and  new  temperature  of  the  ith  piece,  respectively,  and 
the  external  flux  Q  .is  interpolated  to  the  midpoint  of  the  time 
interval.  The  T^4  term  can  be  linearized  in  the  following  way:  expanding 
Ti4  in  a  Taylor's  expansion  about  an  estimated  final  temperature  Ti": 

Ti4  -  T^4  +  4  T"3  (Ti  -  Ti")  (3-37) 

and  substituting  this  expression  for  Ti  in  Ti4  yields  (after  a  further 
approximation) : 
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<*'T 


Ti4  —  [|  <V  +  Tol)  j  4  +  2  (V  +  Toi)J  3  •  (T£  -  Ti')  . 

(3-38) 

Defining  the  term  Tj_'  *  y  (T^'  +  Tq^)  for  convenience,  the  full  linearized 
heat  equation  can  be  written: 


mj  (TP(Tj  *  TP  i  f  j  ,  At\ 

- - - —  *  £/  "J<T  >  «"  -  «.«,3  (t+T)*j 


+2  {Mi}[Ti‘  +  2  Tis  <Ti  -  Tl'>l 

i  * 


Kij 

+  “  <Toi  +  Ti> 


(3-39) 

If  the  terms  in  Equation  3-39  are  grouped  to  isolate  Tj ,  the  following 
equation  results: 


£  c±i  Ti  ° 


(3-40) 


where 


mi  (T/) 


-  2  Mj j  Tj3  -  f  Kj j  for  i  -  j 


-  (2  H±j  T£s  +  \  K±j)  for  i  4  j 


(3-41) 


S' 

r 


*1 


j 


and 


/  At  \  mj  (Tj  "">  Tj  i  r  J 

j  =  Qext,j  (t  +  ~)  Aj  +  At  At  f  “j 

T< 

f  _ a  4  _  Kij  ] 

+  E  [  Mij  (Ti  -  2Ti3  *  Ti0  +  T  ToiJ  • 


B, 


dT 


L°J 


Equation  3-40  is  solved  for  the  Tj  by  matrix  inversion,  and  if  the  Tj  are 
within  one-half  of  one  percent  of  the  estimated  Tj",  the  Tj  are  the  new 
temperatures.  Otherwise,  the  Tj  become  the  estimated  temperatures  Tj" 
and  Equation  3-40  is  evaluated  again  for  T j .  This  process  is  repeated 
until  the  Tj  are  found. 

3.2.2  Estimation  of  At,  T" 

It  is  important  that  the  diagonals  of  the  matrix  c^  be  greater 
than  zero,  so  that  the  time  development  is  stable.  Choosing  the  time 
step.  At,  such  that 


At  < 


®j  (T0j ) 


V  +  1/2  Ki j 


(3-43) 


guarantees  that  the  c^  are  all  positive. 

The  estimated  temperature,  Tj",  is  found  by  solving  the  heat 
equation  with  the  right-hand  side  approximated  by  an  expansion  about  the 
initial  temperature  T0: 

Tj  " 

/  mj  dT  -  At  j  £  [  Mtj  To.“  +  Klj  Toi  +  Q  /  t  +  f\  A 

T«4  (  i  I  V  ' 


+  (4Hjj  Toj3+Kjj)  (Tj'  -  T0j)|  • 


(3-44) 


1 
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Tj is  then  used  to  compute  the  specific  heat  at  Tj',  cpj(Tj'),  and  the 
time  step  is  checked  to  ensure  that 


At  < 


m-jOp 

V  +  Kji  * 


(3-45) 


If  the  time  step  does  not  satisfy  Equation  3-45,  a  smaller  step  is 
chosen,  and  a  new  value  for  is  found  from  Equation  3-44.  This  process 
is  repeated  until  At  is  found  to  satisfy  Equation  3-45.  The  time  between 
trajectory  points  at  which  the  temperatures  are  desired  is  used  as  the 
time  step  if  this  time  is  smaller  than  that  defined  by  Equation  3-45. 

3-3  OPEN  SURFACE  OPTION 

It  is  possible  to  remove  any  station  (more  than  one  station 
may  be  removed)  from  the  replica  or  balloon  shape  to  simulate  an  open 
surface  (hole)  at  that  station.  This  is  accomplished  by  setting  the 
initial  temperature  of  the  station  equal  to  0°R  on  material  property 
card  6.1  (Table  4-3).  The  conductivity  and  radiative  coupling  of  any 
open  station  to  the  rest  of  the  vehicle  is  set  equal  to  zero.  Any 
external  fluxes  entering  the  vehicle  through  the  open  station  are 
neglected,  however,  so  the  percentage  of  open  surface  to  total  vehicle 
surface  should  be  small  to  minimize  the  error  in  determining  the  true 
flux  on  an  interior  surface.  Thus  the  total  number  of  removed  stations 
should  be  small. 

An  average  temperature  is  assigned  to  an  open  surface  for  signa¬ 
ture  calculations  in  BIDIREC.  This  average  temperature  is  representative 
of  the  internal  energy  that  is  passing  through  the  hole  from  the  interior. 

The  average  temperature  of  the  missing  station,  T  ,  is  found  by 

open 

evaluating: 

e  T4  *  Xs  '  g.  T,4  (3-46) 

open  i  ei,open  i 

i 

Because  this  is  an  approximation,  the  internal  emissivities ,  e^,  are 
assumed  to  be  about  the  same,  and  they  cancel  with  the  average  emissivity, 
e.  The  8i>0pen  term  is  the  geometric  factor  for  the  ith  surface  as  seen 
by  the  open  surface. 
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4.  INPUT  SPECIFICATIONS 


A. 1  DISC  UTILIZATION 

MBALL  requires  the  disc  files  used  in  EXOHEAT  for  its  operation. 
Table  4-1  summarizes  the  necessary  tapes  and  their  utilization. 

The  thermophysical  data  (15)  may  be  input  by  the  user.  The 
trajectory  data  (7)  and  earthshine  data  (4)  must  be  input  in  the  BASIC 
option.  Note,  however,  that  the  user  can  supply  trajectory  data  to 
tape  7  by  using  the  trajectory  card  input  option  (Ref.  1). 

TABLE  4-1.  EXOHEAT  (WITH  MBALL)  DISC  FILES 


DEVICE/ 

TAPE 

UTILIZATION 

ORIGIN 

4 

Earthshine  Data 

OSC  Data  Base 

7 

Trajectory  Data 

BALL I S/BASIC 

15 

Thermophysical  Data 

OSC  Data  Base 

16 

Temperature 

MBALL 

23 

Vehicle  Geometry 

RVSNTH/BASIC 

4.2  CARD  INPUT 

The  MBALL  card  input  consists  of  two  parts:  RVSNTH  input  and 
EXOHEAT  input. 

RVSNTH  calculates  the  necessary  geometrical  parameters  for  MBALL 
from  a  modest  user  input  (for  both  replica  and  balloon  shapes).  The  user 
must  be  careful,  however,  to  enter  the  nosecap  polar  angles  [THET(I)]  in 
the  EXOHEAT  input  that  are  calculated  by  RVSNTH  (see  Figure  3-5)  because 
RVSNTH  does  not  pass  them  to  EXOHEAT.  RVSNTH  input  is  displayed  in 
Table  4-2. 
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TABLE  4-2.  RVSNTH  INPUT  SUMMARY 


CARD 

COLUMN 

VARIABLE 

FORMAT 

UNITS 

DESCRIPTION 

FIRST  CARD:  TITLE  (FORMAT  3A4) 

1-12 

(NTI(I), 

I  =  1,3) 

3A4 

RV  title 

SECOND  CARD:  UNIT  DESIGNATOR  (FORMAT  13) 

1-3 

FLAG 

13 

Designates  what  units  of  length 
the  input  is  in: 

_ 

=  1,  ft  =  3,  cm 

=  2,  in.  =  4,  m 

THIRD  CARD:  NOSECAP  (FORMAT  15,  5X,  2F10.4,  15) 

D 

N 

15 

Number  of  stations  on  nosecap 

B 

R1 

FI  0.4 

(FLAG) 

Nosecap  radius  (or  balloon 
radius) 

21-30 

LI 

FI  0.4 

(FLAG) 

Length  of  nosecap*^  (or 
balloon  diameter) 

31-35 

KROL 

’  15 

Number  of  azimuthal  divisions 

FOURTH  CARD:  FRUSTA  CONTROL 

(FORMAT  15,  5X,  F10.4) 

Q 

MN 

15 

Number  of  frusta  (3  maximum) 

B 

Q 

F10.4 

deg 

Cone  angle  of  first  frustum 

FRUSTRA  CARDS:  (ONE  PER  FRUSTUM)  (FORMAT  15,  5X,  2F10.4) 

B 

N 

15 

Number  of  stations  on  frustum 

EES 

Q 

F10. 4 

deg 

Cone  angle 

21-30 

LI 

F10. 4 

(FLAG) 

Length  of  frustum 

BASE  CARD  (FORMAT  F5.0) 

1-5 

AN 

F5.0 

Number  of  stations  on  base 

*See  Figure  3-2  for  replica. 

+If  balloon  shape  is  desired,  LI  =  2  x  R1  and  the  third  card  is  the 
last  RVSNTH  card. 
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The  first  card  contains  the  user's  title  (up  to  12  letters).  The 
second  card  defines  the  units  in  which  the  target  dimensions  are  input 
(these  are  changed  internally  to  feet  for  MBALL) .  The  nosecap  informa¬ 
tion  is  input  on  the  third  card:  N,  the  number  of  stations  refers  to 
the  number  of  longitudinal  stations  along  a  single  running  length  of  the 
nosecap.  The  total  number  of  nosecap  stations  is  KROL  X  N.  For  a  replica 
shape,  the  nosecap  must  fit  tangentially  to  the  first  frustum,  so  that 
LI  is  computed  by  the  user  according  to  Figure  3-2.  If  a  balloon  shape 
is  desired  (i.e.,  sphere),  the  user  should  set  LI  equal  to  twice  Rl. 

RVSNTH  bypasses  the  frustum  calculations  when  LI  ■  2  x  Rl,  so  the 
nosecap  card  is  the  last  RVSNTH  card  when  a  balloon  shape  is  desired. 

For  the  replica  shape,  the  next  card  sets  the  number  of  frusta  and  defines 
the  first  cone  angle,  and  this  card  is  followed  by  a  card  for  each 
frustum.  Again,  the  total  number  of  stations  on  each  frustum  is  N  *  KROL. 
Finally,  the  base  card  contains  the  number  of  longitudinal  stations  on  the 
base  (in  F  format) ,  to  be  multipled  by  KROL  to  get  the  total  number  of 
base  stations. 

EXOHEAT  input  is  shown  in  Table  4-3.  The  first  card  is  the  NAMELIST/ 
SOLAR  data.  This  defines  the  earthshine  mission  map.  The  format  for  this 
card  is  (note  the  leading  blank  denoted  by  the  b) : 

b$SOLAR  ITYP  =  ....  $ 

with  the  variables  separated  by  commas.  A  $  ends  the  namelist. 

The  next  card  set  defines  all  station  normals  by  their  azimuthal 
and  polar  angles.  The  polar  angles,  THET(I) ,  should  agree  with  those 
calculated  by  RVSNTH  (see  Figure  3-5). 

Card  3.1  determines  the  external  flux  averaging  mode  (see 
Table  3-1).  MODE  =  1  allows  the  full  use  of  MBALL’s  capability  for 
azimuthal,  as  well  as  longitudinal,  heat  conduction.  In  MODE  =  2 
(roll  average),  all  azimuthal  stations  experience  the  same  roll  averaged 
external  fluxes.  Azimuthal  conduction  can  still  be  important  in  this 
case,  however,  if  adjacent  stations  have  different  thermal  masses  or 
thermal  properties. 
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TABLE  4-3.  EXOHEAT  (WITH  MBALL)  CARD  INPUT  SUMMARY 


NATURAL  ENVIRONMENT  SPECIFICATION 


NAMELIST/SOLAR 

VARIABLE 

DESCRIPTION 

I  TYPE 

; 

Sun  position  option 

=  1  -  Subsolar  latitude  and  longitude  (SLAT,  SLON) 
input 

=  2  -  Subsolar  point  calculated  from  MONTH,  IDAY, 

HOUR,  where  HOUR  is  GMT  time 

=  3  -  Subsolar  point  calculated  from  MONTH,  IDAY, 

HOUR,  ELG,  where  HOUR  is  local  Sun  time  at 
longitude  ELG 

SLAT 

Sun  latitude  (deg) 

SLON 

Sun  longitude 

MONTH 

Month  of  year  -  1  to  12 

IuAY 

Day  of  month  -  1  to  30 

HOUR 

24  hr  time  (1:15  p.m.  ~  13.25) 

ELG 

Reference  longitude  for  local  time  (deg) 

IPRINT 

Earthshine  map  printing  option 
=  0  -  Do  not  print  earthshine  radiances 
^  0  -  Print  earthshine  map 

I CLD 

Cloud  cover  option 

=  0,1  -  Average  seasonal  geographic  earthshine 

=2  -  Cloudy  radiance  reduction 

=3  -  Clear  radiance  enhancement 

1 

=4  -  Statistical  -  random  geographical  variation 

between  clear  and  cloudy 

=  5*  -  Clear  and  cloudy  sections  are  positioned 
over  the  Earth 

=  6*  -  Same  as  I CLD  =  5,  but  different  albedos  are 
used  for  water,  land,  vegetation,  snow, 
and  ice  on  the  Earth's  surface. 

*See  Table  4-20,  Reference  1. 
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TABLE  4-3  -  Continued 


BODY  SHAPE 

CARD 

PARAMETER 

FORMAT 

DESCRIPTION 

2.1 

KROL 

15 

Numbers  of  azimuthal 
divisions 

2.2  ! 

(ROLL(I),  1=1,  KROL) 

8F10.5 

Azimuthal  angle  (deg) 

2.3 

KTHET 

15 

Number  of  polar  angles 

2.4 

(THET(I) ,  1=1,  KTHET) 

8F10.5 

Polar  angles  (deg)* 

FLUX  AVERAGING  MODE 

3.1 

MODE 

15 

=  1  -  Instantaneous 

=  2  -  Roll  averaged  fluxes 

=  3  -  Spherical  averaged 
fluxes 

BALLOON 

OPTION 

NAMELIST/IBALL 


VARIABLE _ DESCRIPTION 

BALL+  =2  -  Radiative  equilibrium 

=  3  -  Thermal  mass  used 


CONDUCTION  SPECIFICATION 


CARD 

PARAMETER 

FORMAT 

DESCRIPTION 

5.1 

mmm 

T  =  Conduction 

H 

F  =  No  conduction 

COND 

0  =  Thermal  prop,  updated 

II 

with  time 

FLAG 

1  =  Thermal  prop,  not 

■ 

updated 

_ 

*Nosecap  THET(I)  should  agree  with  polar  angles  calculated 
in  RVSNTH  (see  Figure  3-5). 

*BALL  must  equal  2  or  3  for  MBALL  to  be  called. 
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TABLE  4-3  -  Concluded 


MATERIAL  PROPERTY  SPECIFICATION 

(One  card  for  each  station)  The  station  index  for  the  Mth  roll 

angle  and  Nth  polar  angle  is 

I  =  (M-l )  •  KTHET  +  N 

CARD 

PARAMETER 

FORMAT 

DESCRIPTION 

6.1 

JMAT 

(110,  2F10.5) 

Material  identification  code 

word  >200  -  data  taken  from 

OSC  data  base 

<200  data  taken  from  card 

sets  composed  of  card  types 

7.1  and  7.2 

THK 

Material  thickness  (ft) 

TOLD 

Initial  temperature*  (°R) 

If  some 

of  JMATS  are 

less  than  200  material  property  data  will  be 

input  on 

these  cards 

7.1 

NCP 

(15,  7F10.2) 

Number  of  temperatures  and 

number  of  7.2  cards 

DENS 

Density  (lb/ft2) 

EO 

Outside  emissivity 

El 

Inside  emissivity 

AL 

Absorptance 

m 

TEMP 

(3E12.5) 

Temperature  (°R) 

CP 

Specific  heat  (Btu/lb/°R) 

TCON 

Thermal  conductivity 

■ 

_ 

(Btu/ft/°R/sec) 

♦TOLD  =  0  for  an  open  surface 
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The  next  card  is  NAMELIST/IBALL  input.  If  BALL  *  2,  a  zero 
thermal  mass  is  assumed  for  all  of  the  stations,  and  Equation  3-34  is 
solved  for  the  temperatures.  BALL  ■  3  uses  the  complete  thermal  mass 
solution  of  Section  3.  BALL  must  equal  two  or  three  for  MBALL  to  be 
called.  The  format  for  this  namelist  is: 

b$IBALL  BALL  -  2 (or  3)  $ 

Card  5.1  chooses  the  two  options:  conduction  on/off,  and  thermal 
properties  updated/not  updated  with  time. 

Card  set  6  defines  the  material,  thickness,  and  initial  tempera¬ 
ture  of  each  station  (one  card  per  station) .  An  open  station  is  flagged 
by  setting  its  initial  temperature  equal  to  0.  If  JMAT  is  greater  than 
200,  the  thermophysical  properties  are  taken  from  the  OSC  data  base. 

Table  4-4  contains  the  OSC  thermophysical  property  code  words.  If  JMAT 
is  less  than  or  equal  to  200,  the  thermophysical  properties  must,  follow 
card  set  6.  If  a  station  is  built  up  of  layers  of  different  materials, 
an  appropriate  average  of  their  thermophysical  properties  over  the 
thickness  of  the  skin  should  be  input  in  card  set  7  (see  Equations 
3-31  through  3-33) .  Card  set  7  is  read  whenever  a  JMAT  less  than  or 
equal  to  200  is  encountered  that  is  different  from  the  previous  station's 
JMAT. 

To  model  a  phase  change,  the  specific  heat  can  be  thought  of  as 
a  spiked  function  of  the  temperature,  as  is  shown  in  Figure  4-l(a).  The 
enthalpy,  or  energy  content  per  unit  mass  of  the  station  material, 
is  the  area  under  curve  (a) ,  and  is  shown  in  (b) .  The  area  of  the 
shaded  spike  of  curve  (a)  corresponds  to  the  heat  of  fusion  or  heat 
of  vaporization  of  a  unit  mass  of  material.  In  the  case  of  water,  this 
area  would  correspond  to  about  144  Btu/lb  for  the  heat  of  fusion. 

Table  4-5  is  a  tabulation  of  the  curve  of  Figure  4-l(a).  Note  that 
the  width  of  the  spike  was  chosen  to  be  1°R,  but  a  smaller  width  may 
be  chosen  if  the  Cp  is  increased  so  that  the  area  under  the  spike 
remains  144  Btu/lb. 


ENTHALPY  (Btu/lb) 


TABLE  4-4.  OSC  THERMOPHYSICAL  PROPERTY  CODE  WORDS 


CODE 

WORD 

MATERIAL 

■■ 

Carbon  phenolic 

Graphite 

Kfl 

Silica  phenolic,  asbestos  phenolic 

204 

Fused  silica 

205 

Teflon 

206 

Porous  stainless  steel 

207 

Aluminun 

208 

Beryllium 

TABLE  4-5.  MODELING  THE  HEAT  OF  FUSION  OF  WATER 


TEMPERATURE  (°R) 

cp(Btu/lb/°R) 

480.00 

0.5  (approx.) 

491.49 

0.5  (approx.) 

491.50 

144.5 

492.50 

145.0 

492.51 

1.0  (approx.) 

500.00 

1.0  (approx.) 

4.3  ACCESSING  MBALL  IN  THE  BASIC  OPTION 

MBALL  Is  accessed  in  the  BASIC  option  (note  that  Reference  1 
and  3  do  not  contain  MBALL)  if  the  BASIC  parameters  HEATRV  *  EXOHEAT 
and  TARGSYN  ■  YES  (this  calls  RVSNTH) .  References  1  and  3  state  that 
RVSNTH  should  not  be  used  with  EXOHEAT,  but  this  has  been  changed  with 
the  addition  of  MBALL.  To  compute  signatures  from  the  MBALL  temperatures, 
input  is  necessary  for  SELECT  and  BIDIREC,  because  the  modified  RVSNTH 
program  does  not  compute  the  necessary  parameters  to  these  programs. 
Reference  1  describes  SELECT  and  BIDIREC  input. 
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5.  EXAMPLES 


5.1  WATER-JACKETED  BALLOON 

The  first  example  (Figure  5-1)  demonstrates  the  modeling  of  a 
water- jacketed  balloon  1  m  in  radius  consisting  of  an  outer  and  inner 
skin  of  1/32-in.  teflon,  supporting  a  1/8-in.  layer  of  water.  It  is 
necessary  to  average  the  density,  specific  heats,  and  thermal  conductivities 
of  the  water  and  teflon  over  the  layers  of  the  balloon.  The  balloon's 
surface  is  divided  into  36  stations  of  equal  area.  Table  5-1  shows  the 
values  chosen  to  represent  the  specific  heat  of  water  to  model  the  phase 
change  at  492  °R,  and  the  specific  heat  of  teflon.  Also  shown  are  the 
respective  thermal  conductivities.  The  averages  used  in  the  inputs  were 
determined  from  Equations  3-31  through  3-33. 

Trajectory  cards  have  been  input  to  place  the  balloon  over  the 
north  pole  with  the  axis  of  the  balloon  parallel  to  the  axis  of  the  Earth. 
The  Sun  is  at  0°  longitude  and  0°  latitude.  The  initial  temperature 
of  the  entire  balloon  is  500  °R.  Even  though  SELECT  and  BIDIREC  input 
is  present,  the  output  from  these  programs  is  suppressed. 

The  MBALL  output  (Figure  5-2) ,  following  the  average  fluxes  from 
the  Sun,  Earth,  and  molecular  heating  from  EXOHEAT,  is  as  follows  for 
each  station  (units  are  in  feet):  the  station  area,  perpendicular 
and  parallel  components  of  the  station  normal  with  respect  to  the 
vehicle  axis,  distance  from  the  axis  (r) ,  distance  along  the  axis  (z) , 
running  length,  and  thermal  mass  (Btu/°R).  The  terms  RAD  and  COND 
give  an  estimate  of  the  initial  radiation  and  conduction  flow,  respectively, 
involving  the  particular  station.  TMTP  is  the  estimated  time  step  At 
given  by  dividing  the  thermal  mass  TMASS  by  the  sum  of  RAD  and  COND  (this 
is  one-half  the  At  of  Equation  3-43).  ALF,  EOUT,  and  EIN  are  the  solar 
absorptivity  and  outside  and  inside  emissivities ,  respectively.  The 
temperatures  (°R)  are  then  output  for  each  station. 
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FIGURE  5-1.  BALLOON  EXAMPLE  INPUT 
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TABLE  5-1.  p,  cp,  k  FOR  WATER  AND  TEFLON 


p 

(lb/ft3) 

CP 

Btu/lb/°R) 

k 

Btu/ft/°R/sec) 

T(°R) 

62.4 

0.5 

0.323 

480.00 

62.4 

0.5 

0.329 

491.49 

62.4 

144.5 

0.329 

491.50 

62.4 
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492.50 
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0.334 
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0.414  x  io-4 
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0.414  x  io’4 

491.49 

135.0 
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0.414  x  10"4 

491.50 
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0.414  x  IO"4 
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135.0 
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0.414  x  IO-4 

492.51 

135.0 
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5.2 


REPLICA 


This  example  consists  of  a  biconic  shape  with  a  tangentially 
fitting  spherical  nosecap  that  is  divided  into  24  stations,  eight 
longitudinal  by  three  azimuthal.  The  nose  is  made  of  silicon  phenolic 
that  tapers  from  0.3  in.  for  the  first  two  nose  stations  (along  a 
single  longitudinal  ray)  to  0.2  in.  for  the  third  nose  station.  The  first 
frustum  has  two  stations,  each  consisting  of  a  0.1- in.  aluminum  structure 
covered  by  a  0.1-in.  thickness  of  silicon  phenolic.  The  second  frustum 
contains  two  stations  with  the  a  0.1- in. -thick  aluminum  structure  covered 
by  0.05  in.  of  silicon  phenolic.  The  base  has  one  station  of  material 
structured  similar  to  that  of  the  first  frustum.  The  dimensions 
of  the  replica  are  shown  in  Figure  5-3.  Note  that  whenever  the  material 
identification  code  word  JMAT  of  card  6.1  changes,  card  set  7  must  be 
input.  The  nose  is  completely  silicon  phenolic  (JMAT  «  203),  and  the 
numbers  JMAT  «  180  and  JMAT  =  190  have  been  arbitrarily  assigned  to  the 
aluminum  and  silicon  phenolic  combinations.  Equations  3-31  through  3-33 
were  used  to  compute  the  average  values  of  p,  cp,  and  k  that  are  input. 
Figure  5-4  and  5-5  show  the  input  and  output,  respectively,  for  this 
example. 
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NOTE:  ALL  DIMENSIONS  IN  cm 


FIGURE  5-3.  REPLICA  GEOMETRY 
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FIGURE  5-4.  REPLICA  EXAMPLF  INPUT 
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FIGURE  5-4  -  Continued 


exoHE*r  hill  c*lCul*T€  thermal  not* 


iNNNrgotNSi^NAJMNN 


'tJW4^WUiuJL.k 


i  uu  u  u  jwWU'wwu  uj u> 


0»1o«#*dQ«aaQ^aBttit  *40  00 


*4T4«4*4f1>4(VNNNM(^Mrvj(VNMNNNNNIW(V 

U/WWUMU/'iiWWUllUltiU.'WUJItMlIWUJItlWWWW 


%D 


CiC^ivocS^NClCifi  ?  N  a 
iwM^miCifitkCir  a  tun  4 


UbiUtU-'UJtMb'WUlUtbUtUlU*  1IWWUI1.  ujtMUJku.' 
a  on  o  o  a  r'  o  o  »•  o  o  csnuonnca  ere 
Dooooonuo j  o • ^ oo  a*  n  u  bo  430 
'jraae^oooooQooodc  9  7 cq  ■  ac 

(tC 

irtn«4trtwtfi<i»n«4Nt*)«iruv<oNc09  9*1  m  tn  n 
if'  >fl  K  »*  « 


oaaooArnaoooooas*  j  c.  o  a  e  a  • 

**♦♦♦♦♦♦♦♦♦♦*+*♦♦******♦ 

UtU-tilUIWUlhJU'MJWUiW  W^uiUJUlWbJUJUJ  MIW<al 

•>ao00Oaaaooao«UJofrauaaor 
naaoattcoott*  9  o  s*  >a  o  c  a  r  aOu 

C.  30ULltfCh*UUL«as4s«Ulrf0>*C  bo  JILtd 
O  93  0OO  (MOC  0  0NCt«<£  J«iOCifi4  M«C 
jr  NO»iC4c^^(v«/u'UMfiN«(re  o-^Mnn 


vt^V4«4V4<4»tjNNNN(LlNN(VOtlVJNNN(yj(yM\tN 

a  J>aoCriwadaaa.jwaaab.t.  a  a  an 
#♦♦♦♦**♦♦*♦*♦*♦♦«*  +  ♦*+*♦ 
UJ  Ui  uJ  Ui  4*1  Ul  Ui  Mi  tei  IUH<  Ui  II)  b1  ui  III  Itl  UJ  b>  Id  UJ  U  U  «•» 

aooouoevovooanaaoaau  a  oji 


a  5  a  &  a  *  t-  9  a  ac  0  a  t  a  3r  c  n  r  w  r  o  a 
+  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  «-***♦ 
UiU.UWUIWWk  WitiuMi  UI  U  b.  li  (it  U1  tb  w  U4  Ul  h/  1> 
u  aaocaoooot  aaouaaau  a  r.  ^  o  u 

Qi  ajooaooaaoouooaua  oar 

•1i)0«O9«»iAINaC  iT  -f  C  C  if  J(VJ  U«  4> 
rw  f»m  if  fNOi.HM'lJ  0  UMf  S  C  c  (T  w  *4  <V  (M  P) 

IT  1010  N  C(f  «<*(««  «4«4«4  ««M(V  (VN<V 


«4«4«4«4*4«4«4(V(V<WriJCW<Mr\j04fVJ(M<^rv<Vf\J  .aOICNa 
aaaaoff wo09 1 >000040001 000 
*♦♦♦♦♦*♦♦*♦♦♦♦♦♦♦♦♦♦♦  ♦♦♦ 
UlUJUiUltaJUJUiUiUlWbIWWUiUlUJUlfeJUi  -JUl  Wi  UJ  UI 

o An JooorAoereaoaooaAw ooo 

OOOAOOOOOACtOOOOVOOf  10  3  0  0 

6*ttOJBOM(t«*»0MM)*<»KlM<l^  *r  K.  IIS 

HS'Kifn*«7a«4eunn>00K.Ko<raH0N0 


•mO«(«(«(ftNllNNNAMWNNMN(MNMMMN 

0  Ui|i/Ufli|Ut|b  ultlllil  ItfkituutWItlUIUIllW  uilu  dlWUI 

DBOOO0O0OOOO  >ooo«ooao  ooo 

«  «t/OO0»»,'04  N  »t  C  J  NCd  if  4  NU  C  if  4 

s  uc4«K9«o«4AnA#00Mh«0'O*4««on 

-J  r  tft  <*  *•  •  *  r  •*  *4  «•  «4  ««  «4  «4  •«  ««  v«  «4  *4  MM  *I«U  V* 


*HI6  PAOl  Js  Eiv  T  QWjI27  rjy  -Tmvrr- 


FIGURE  5-5.  REPLICA  EXAMPLE  OUTPUT 
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FIGURE  5-5  -  Concluded 


6.  HINTS  AND  DIAGNOSTICS 


To  call  MBALL,  it  is  necessary  that:  HEATRV  =  EXOHEAT 

TARGSYN  -  YES 

BALL  =  2  or  3  (in  NAMELIST/ 
IBALL) 


If  the  removal  of  a  piece  (creation  of  a  hole)  is  desired: 

a.  The  initial  temperature  of  the  piece  is  input  as  0. 

b.  The  external  flux  entering  the  interior  is  neglected,  so 
the  number  of  stations  removed  should  be  small. 

For  a  skin  that  consists  of  layers  of  different  materials,  the 
appropriate  averages  of  the  material  properties  should  be  input 
(see  Equations  3-31  through  3-33)  in  card  set  7  (Table  4-3) . 

If,  for  any  reason,  it  is  necessary  to  model  a  station  with 
an  outside  emissivity  equal  to  zero  (this  piece  thus  has 
no  radiative  communication  with  the  external  environment)  the 
outside  emissivity  on  card  7.1  (Table  4-3)  should  be  set  to  a 
small  number,  but  not  zero.  This  is  because  the  external 
emissivity  is  set  equal  to  zero  internally  for  an  open  station 
(when  TOLD  =  0)  and  is  used  as  the  flag  for  that  open  station 
in  the  calculations,  so  that  any  outside  emissivity  equal  to 
zero  would  be  interpreted  as  an  open  station. 
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